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Abstract 


^  In  a  packet  switched  network,  several  types  of  transmission  may  connect  a  pair  of 
nodes  of  the  network.  The  delay  characteristics  of  each  type  of  transmission  may  be 
significantly  different  and  so  a  single  queue  to  the  node_pair  operating  under  a  first- 
come,  first-served  rule  may  not  be  best.  Mn-this  paper^we  develop^both  exact  and  ap¬ 
proximate  mathematical  models  that  determine  various  system-level  performance 
measures  for  a  (K,N)  scheme.  Under  this  scheme  arriving  packets  are  placed  in  the 
queue  for  one  type  of  transmission  until  the  size  reaches  K,  whereupon  arriving  pack¬ 
ets  are  diverted  to  the  other  transmission  system.  When  the  number  of  packets  at  the 
first  system  drops  to  N(<K),  arriving  packets  are  again  placed  in  the  first  queue.  / 


I.  Introduction 


In  a  packet  switched  communication  network,  several  transmission  paths  may  exist 
between  two  switches  in  the  network  (see  [1]).  The  delay  characteristics  of  each  path 
may  be  significantly  different;  for  instance,  one  may  be  a  satellite  link  (with  its  inher¬ 
ent  propagation  delay)  and  the  other  a  56kb/s  terrestial  link.  Thus  if  packets  from  the 
same  message  are  processed  at  the  two  links  on  a  first-come,  first-served  basis,  they 
will  be  received  at  the  next  node  out  of  order.  In  order  to  alleviate  this  problem,  a 
(K,N)  scheme  can  be  formulated  for  the  operation  of  the  system. 

We  consider  the  situation  where  there  are  two  types  of  transmission  and  arriving 
packets  are  placed  in  the  queue  for  the  type  1  transmission  system  until  the  number  of 
packets  in  that  system  is  K,  whereupon  packets  are  sent  to  the  second  transmission 
system.  Packets  continue  to  be  sent  to  the  second  queueing  system  until  the  number  of 
packets  at  the  first  queueing  system  drops  to  N(<K),  at  which  time  arriving  packets  are 
again  placed  in  the  first  system.  The  following  figure  depicts  this  scheme,  which  is 
known  as  the  of  heterogeneous,  multi-trunking  queueing  system. 


Figure  1.  Heterogeneous  Multi-Trunking  Queueing  System 

We  are  considering  a  queueing  system  composed  of  two  queues  and  a  single  arri¬ 
val  process  of  customers  (packets).  Arriving  customers  are  placed  into  one  or  the  other 
queue  based  on  the  number  of  customers  at  the  first  queue.  Once  assigned  to  a  queue 
they  remain  there. 

Queueing  systems  similar  to  this  one  have  been  considered  before.  In  [2]  Singh 
examined  the  single  queue,  first-come  first-served  version  of  the  system.  That  is,  cus¬ 
tomers  are  placed  in  service  when  a  server  is  free  based  on  their  order  of  arrival.  An¬ 
other  version  of  this  system  is  the  shortest  queue  problem  where  arriving  customers 
are  placed  in  the  queueing  system  with  the  smallest  number  of  customers  (see  [3]  and 
(4)  and  the  papers  referenced  therein). 

In  Section  II  we  give  a  formal  definition  of  the  problem,  our  assumptions  and  a 
complete  mathematical  analysis  of  the  system.  The  analysis  required  finding  roots  of 
an  equation  and  the  solution  of  a  set  of  linear  equations.  In  light  of  this,  we  developed 


an  approximation  to  the  system-level  performance  that  does  not  require  these  poten¬ 
tially  complicated  numerical  problems.  Section  lH  documents  the  approximate  model. 
Numerical  comparisons  between  the  exact  and  approximate  models,  as  well  as  other 
numerical  examples,  are  given  in  Section  IV.  Finally,  Section  V  gives  a  summary  of 
the  paper. 


II.  Mathematical  Performance  Model 


We  are  considering  the  system  shown  in  Figure  1.  We  assume  customers  (packets) 
arrive  in  accordance  with  a  Poisson  process  with  rate  A..  Let  Qp  i=l,2,  be  the  steady- 
state  number  of  customers  in  queueing  system  i.  If  Qi=K,  then  arriving  customers  are 
placed  in  queueing  system  2.  This  continues  until  Qi=N  (N<K),  whereupon  arriving 
customers  are  again  placed  in  queueing  system  1.  Once  placed  in  a  system,  a  customer 
remains  there.  We  assume  that  the  length  of  time  to  serve  a  customer  at  server  i  is 
exponentially  distributed  with  mean  1/p.j,  i=l,2,  and  the  service-time  distributions  and 
arrival  process  are  all  independent  of  each  other. 

Under  this  rule  the  number  of  customers  in  system  1  is  at  most  K,  whereas  there 
could  be  an  infinite  number  at  system  2.  The  value  of  N  that  is  currently  being  consid¬ 
ered  in  packet  switched  networks  is  N=K-1,  so  we  first  present  an  analysis  of  that 
case. 

Case  II. a:  The  (K,K-1)  System 

In  this  system,  when  the  number  of  customers  in  system  1  reaches  K,  arriving 
customers  are  sent  to  system  2  until  we  have  the  first  departure  from  system  1.  Let  us 
begin  the  analysis  by  defining 

Pij  =  Pr{QH,  Q2=j}  (i=0,l . K;  j=0,l . ).  (1) 


Then  the  behavior  at  system  1  is  given  by  the  results  of  an  M/M/l/K  queueing  system. 
Let  Pj  =  xyjij,  i=l,2;  then  from  Gross  and  Harris  [5],  we  have  for  i=0,l,...,K, 
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The  steady-state  equations  for  p;.  are 
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For  i=0,l,.--,K,  and  |z|^l  define 
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pi (z)  =  £  PjjZj  . 

j=o  1 

Then  we  have  the  following  matrix  equation  for  the  Pj(z): 
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Systems  of  equations  like  these  appear  frequently  in  these  types  of  queueing  problems 
and  the  solution  technique  is  standard;  for  instance,  see  [6],  [7],  and  [8]. 

Define  the  following  recursions: 


and 


E-i(z)  =  1 
E0(z)  =  a(z) 

Ej(z)  =  b(z)Ei_1(z)-X.p1z2Ei_2(z) 
Ek(z)  =  c(z)Ek_1(z)  -kp1z2EK_2(z) 


(7) 


(i  <K-D 


F_i(z)  =  1 

E0(z)  =  c(z)  (8) 

Fj  (z)  =  b(z)Fi_1(z)  p1z2Fi_2  (z)  (i  <K-1) 

FK(z)  =  a(z)FK_1(z)  -\p1z2FK_2(z). 

If  A(z)  is  the  (K+l)x(K+l)  matrix  multiplying  the  column  vector  of  (Pj(z)}  on  the  left 
hand  side  of  Equation  (5),  then 


det[A(z)]  =  Ek(z)  =  Fk(z). 


(9) 


Now  let  Aj(z)  be  the  matrix  A(z)  with  the  right-hand  side  of  Equation  (5)  in  the  ith 
column,  then,  for  i=0,l,...,K  we  have 


P,(z)  =  det'A>(z)l  , 
det[A(z)] 


(10) 


where 


i-l 


det[Aj  (z)]  =  P-2  (z-1)  S  (\z)'"k  Pk0Ek_i(z)FK-i-i(z) 


(ii) 


with  the  first  sum  set  to  zero  when  i=0.  Thus  the  problem  is  solved  once  we  find 
p.Q,  i=0,..., K. 

We  need  to  develop  K+l  equations  in  the  unknowns  {p.Q  i=0,l,...,K}.  The  first  can 
be  obtained  from  the  normalizing  condition  (  2P.(l)=l).  It  is  straightforward  to 

v  i=0  1  7 

show,  using  this  condition,  that 


i=0 


(12) 


i-P“ 


K+l 


and  that  the  right-hand  side  of  Equation  (12)  must  be  positive  for  Q2to  be  a  well  de¬ 
fined  random  variable.  Since  the  quantity  p*(l-p  )/(l-p*+1)  is  the  portion  of  time  cus¬ 
tomers  are  arriving  at  system  2,  p2  times  this  quantity  is  the  effective  load  on  system 
2,  and  this  has  to  be  less  than  1  in  order  for  p..,  to  be  well  defined.  Thus  the  require¬ 
ment  that  the  right-hand  side  of  equation  (12)  be  positive  makes  intuitive  sense.  There 
are  no  existence  conditions  for  Qlt  unless  K=oo. 

It  is  easy  to  show  that  det[A(l)]ssl  and  so  z=l  is  a  root  of  the  denominator  of 
Equation  (10).  Using  a  similar  analysis  as  in  [6]  or  [8]  one  can  show  that  there  are  K 
distinct  roots,  £.,  i=l,2,...,K,  of  det[A(z)]=0  in  (0,1).  Since  the  numerator  of  Equation 
(10)  must  also  vanish  at  these  points,  we  have  K  additional  equations  given  by 


i  («l)K'kPkoEk-l(«i)-°  0-1.2 . K).  (13) 

k=0 

Thus  Equations  (12)  and  (13)  give  us  the  required  K+l  equations  in  the  {p.  Q} . 
Summing  Equation  (5)  we  also  have 
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where  PK(z)  is  given  by  Equations  (10)  and  (11).  From  Equation  (14)  we  have 
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Using  Equations  (2)  and  (14)  one  can  compute  the  expected  values  of  Qi  and  Q2 
and  hence  the  expected  value  of  the  number  of  customers,  OQ1+Q2,  in  the  system. 
The  average  waiting  time,  E{W},  is  then  given  via  Little’s  Formula,  E{W}  =  X_1E{Q}. 
For  the  case  when  K=1  and  N=0,  we  have 


(\  +  M.2)z-P2  ~^\z 

P0(z> 

Poo 

=  P2(z-1) 

~Kz  -X  z  +(X.  +  |x1  +  p2)z -p2 

_Pl  ( zl 

_P  10_ 

(16) 


det[A(z)]  =  (z— 1)  (~\2z2  +  i2)z-ji2) 
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Even  in  this  simple  case  no  nice  expression  exists  for  E{Q2},  so  one  is  forced  to  dif¬ 
ferentiate  P(z). 

Case  II. b:  The  (K,N)  System 

For  the  case  of  a  general  N,  we  have  to  redefine  p^  to  take  into  account  which 

system  is  receiving  arriving  customers.  When  N=K-1  the  value  of  Qi  determines  which 
system  receives  arriving  customers.  Begin  here  by  letting  I  be  an  indicator  variable 
such  that 

10  when  arriving  customers  go  to  system  2 
1  when  arriving  customers  go  to  system  1 

and  define 

P.!  =  Pr{Qi=i,  Q2=j,  1=1}- 

ij 

Now  if  1=0  then  p°=0  for  i=0,l,...,N  and  all  j.  If  1=1  then  p*  =0  for  all  j.  As  in  case 
ij  Kj 

Il.a  one  can  write  down  the  steady-state  equations  in  p1  and  the  resulting  generating 
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and  a(z),  b(z),  and  c(z)  are  defined  as  before.  We  have  annotated  in  which  of  the 
rows  and  columns  of  A(z)  structural  changes  occurred.  Again  we  have 
(P°(z)  =  0  for  i=0,l,...,N) 

P°(z)  =  M-2(z-D  det[Ai(z)]  (i=N+l,...,K);  (24) 

1  det[A(z)] 

where  A,(z)  is  the  matrix  A(z)  with  P0  in  the  K+i-N  column.  We  start  numbering  the 

columns  of  A(z)  with  0.  For  1=1  we  have  an  analogous  form  for  p|(z)  for  i=0,l,...,K-l 

as  that  given  for  P^(z),  except  that  det[Aj(z)]  is  the  matrix  A(z)  with  P0  in  the  \th  col¬ 
umn.  There  is  no  simple  recursion  we  can  give  for  det[A(z)]  or  det[Aj(z)]  as  was 
done  in  the  case  of  N=K-1.  But  the  analysis  is  still  the  same;  namely,  one  has  to  find 

the  2K-N-1  roots  of  det[A(z)]  in  (0,1),  and  set  the  numerator  of  one  p'(z)  to  zero  at 
these  roots  to  find  2K-N-1  equations  for  p^Q,  i=N+l,...K  and  p.1^  i=0,l,...,K-l. 

The  other  equation  is  obtained  from  the  normalization  condition  and  we  have 

Z  p°  +  Z  p1  =  l-p  Pr{I=0} .  (25) 

i=N+l  10  i=0  10  2 


Thus  we  need  only  find  an  expression  for  Pr{I=0}  to  complete  the  solution. 
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In  order  to  develop  an  expression  for  Pr{I=0}  we  use  first-passage  results  from  the 
standard  M/M/1  queueing  system  (see  [9]).  Let  t-  be  the  average  first-passage  time  in 

system  1  from  state  i  to  state  j;  then 


Pr{I=0}  = 


tkn  +  tnk 


In  system  1  the  expected  time  to  go  from  state  K  (Qi=K)  to  N  (Qi=N)  is  the  time  for 
K-N  departures,  so  that 

Tkn  =  (K-N)/pi.  (27) 

Now  to  go  from  state  N  to  state  K  for  the  first  time  in  system  1  is  the  same  as  in  a 
standard  M/M/1  queue,  and  is  given  in  [9]  as 
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Using  all  of  these  results  we  thus  have 
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We  may  also  study  the  behavior  at  system  1  in  isolation.  Using  z=l  in  the  equa¬ 
tions  for  pj(z)  one  gets  the  following  equations  for  p{  =  Pr{Qi=i,  1=1}: 

m>o-h,p!  =  ° 
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(30) 
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Pn+i=Pn+2  =--  =  Pk  =  PiPk-i 


The  solution  to  these  equations  is 
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from  which  we  have 
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For  the  case  when  pi=l  we  have 
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III.  Approximate  Behavior  at  System  2 
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In  Section  n  we  gave  a  complete  analysis  of  the  behavior  of  system  1,  with  results 
contained  in  Equations  (32)-(36).  The  behavior  of  system  2  rests  on  finding  the  roots 
of  an  equation  and  then  solving  a  system  of  equations.  The  results  for  system  1  are 
computationally  simple  enough;  it  would  be  nice  to  have  similar  simple  expressions  for 
system  2.  Its  behavior  is  similar  to  a  standard  M/M/1  queueing  system  where  the  arri¬ 
val  process  of  customers  alternates  between  on  and  off  times.  Kuczura  (10]  presented 
the  idea  of  an  interrupted  Poisson  process  (IPP)  as  a  way  of  analyzing  overflow 
streams  in  loss  systems.  Heffes  [11]  extended  this  work  to  a  GI/M/1  queueing  system 
where  the  arrival  process  is  an  interrupted  Poisson  process.  His  results  will  be  used 
here  to  approximate  the  behavior  at  system  2.  For  an  IPP  it  is  assumed  that  the  on 
and  off  time  of  the  underlying  Poisson  process  is  described  by  independent  exponen¬ 
tially  distributed  random  variables.  In  general,  this  would  not  be  true  at  system  2,  ex¬ 
cept  in  the  case  of  K*l,  N=0.  But  the  results  of  Heffes  should  provide  an  excellent  ap¬ 
proximation  to  the  behavior  of  Q2. 

More  formally,  for  system  2  we  have  a  basic  Poisson  arrival  process  with  rate  k. 
The  input  is  turned  on  (off)  for  a  period  of  time  that  is  exponentially  distributed  with 

mean  y-1  (w-1).  Following  Heffes  [11]  the  interarrival  distribution  at  system  2  of  cus¬ 
tomers  is  the  mixture 

A(t)  =  ki  (l-e~rit)  +  k2(l-e_r2t)  ,  (37) 


where  for  our  system  we  have 


y  =  Hi  /(K-N), 


A(l-P, ) 


K-l  l-pn+1  ’ 

y  _  Mi 


n  n 

n=N  pi 


rj  =  1  {  \+w  +7  +  \/(\+u+y)2  -  4^o>}  , 


r2  =  1  {  +y-V(\+w+7)2  -  4\w}  , 


i»;v. 


k2  *  1  -  k, 


(42) 


A. -r. 


ri-r2 


Using  standard  results  from  the  GI/M/1  queueing  system  (see  [5]),  we  have  the 
approximation  to  the  number  of  customers  at  system  2  given  by 


Ea{Q2}  = 


P2  Pr{I=0) 
1-0 


where  Pr{I=0}  is  given  by  Equation  (29)  and 


(43) 


0  = 


_~T]2  W  ^2  ~  ^1^3 

2t){ 


with 


(44) 


Til  =  ~V\, 

t|2  -  n  42  +  n  P2  +  and 
TJ3  =  -(ki  n  P2  +  k2  r2  P2  +  n  r?)- 

We  have  just  written  down  the  expected  value  of  Q2;  but  using  the  results  from 
the  GI/M/1  queue  one  can  give  an  expression  for  an  approximation  to  the  probability 
distribution  of  Q2  as 


Pr{02  =i}  = 


1  -  p2  Pr{I=0) 
p2  Pr{I=0>  (1- 0)9  j'1 


0=0) 


0=1,2,...). 


(45) 


One  can  use  this  approximation  along  with  the  exact  results  for  system  1  (see 
Equation  (35))  to  give  an  easily  computed  approximation  for  the  total  number  of  cus¬ 
tomers  in  both  systems.  The  average  waiting  time  can  then  be  found  via  Little’s  For¬ 
mula.  In  the  next  section  we  compare  the  exact  results  with  approximations. 
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IV.  Numerical  Examples 

In  this  section  we  present  some  numerical  examples  using  the  results  that  were 
developed  in  Sections  II  and  HI.  We  investigate  three  areas: 

a.  Comparison  of  exact  and  approximate  models; 

b.  Characterization  of  optimal  K  in  the  (K,  K-l)  system  that  minimizes  expected 
number  of  customers  in  the  system;  and 

c.  Comparison  of  multi-trunking  with  first-come,  first-served  schemes. 

Table  1  presents  some  comparisons  of  the  exact  results  for  the  case  K=2,  N=l, 
and  the  IPP  approximation  of  Section  HI.  In  the  IPP  case  we  are  approximating  the  on 
and  off  periods  of  the  arrival  process  to  system  2  by  exponential  distributions  with  ap¬ 
propriate  means.  Since  we  can  give  exact  values  for  the  mean,  the  only  facet  of  the 
approximation  comes  from  the  form  of  the  underlying  distribution.  When  K=1  and  N=0, 

the  on  and  off  times  are  exponentially  distributed  with  means  p.”1  and  A.-1  and  the 

approximate  and  exact  results  do  agree.  For  the  case  of  N=K-1,  the  on  time  is  again 
exponentially  distributed  with  mean  1/p.i  and  so  only  the  off  time  distribution  is  being 
approximated.  Thus  for  the  major  cases  of  interest  one  would  expect  the  approxima¬ 
tion  to  be  good.  This  is  evident  in  Table  1. 


^1  =  ^  =  1 

=  5,  n,  =  1 

=  10,  ^  =  1 

n 

E{Q} 

n 

E{Q} 

D 

E{Q} 

.3 

.3682 

.3680 

2 

.8073 

.7944 

2 

.3082 

.3068 

.6 

.8320 

.8287 

2.5 

1.352 

1.316 

3 

.6556 

.6461 

.9 

1.432 

1.415 

3 

2.467 

2.367 

4 

1.414 

1.372 

1.2 

2.402 

2.346 

3.5 

6.213 

5.863 

5 

4.168 

3.974 

TABLE  1.  Comparisons  of  Exact  and  Approximate  Models. 


In  Table  1  we  have  compared  the  results  of  the  approximation  with  the  exact  re¬ 
sults  for  the  case  K=2,  N=l.  The  exact  results  are  given  as  the  upper  entry  in  each 
box.  For  increasing  arrival  rates  three  comparisons  were  made;  namely, 

M-i/p-2  =  1.  5,  and  10.  For  each  case  considered  we  see  that:  the  approximation  under¬ 
estimated  the  exact  result;  as  the  arrival  rate  increased  the  approximation  got  worse; 
and  the  ratio  of  the  service  rates  was  not  a  factor  in  the  goodness  of  the  approxima¬ 
tion. 

Since  the  approximation  is  underestimating  the  exact  results,  one  would  expect 
that  the  off  time  of  the  arrival  process  at  system  2  would  be  better  approximated  by  a 
distribution  whose  coefficient  of  variation  is  greater  than  one.  We  know  of  no  simple 
results  for  the  IPP  with  non-exponential  on/off  times.  The  best  the  approximation  did 
was  for  the  case  X  =.3,  pi  =  p.2  =  1,  with  a  relative  percent  error  of  .054%;  and  the 
worst  case  was  a  5.6%  error  for  X.  =  3.5,  p.i  =  5,  and  p.2  =  1-  One  would  expect  this 
type  of  behavior  because  as  X.  increases,  the  load  on  system  2  increases  and  the  ap¬ 
proximation  is  being  used  to  a  greater  extent.  But  even  in  the  case  of  X.  =  3.5,  pi  =  5, 
and  P2  =  1,  the  accuracy  is  quite  good.  The  fact  that  there  were  not  any  great  differ¬ 
ences  in  the  accuracy  of  the  approximation  as  pi/p.2  ranges  from  1  to  10  is  quite  a  sur¬ 
prise  and  perhaps  follows  from  the  fact  that  the  mean  off  time  (see  Equation  (39))  de¬ 
pends  only  on  X.  and  pi  and  not  on  pi  and  p.2-  All  in  all,  we  are  quite  satisfied  with 
how  well  the  EPP  approximation  does. 

When  cost  is  not  a  consideration,  one  might  guess  that  the  best  value  of  K  in  a 
(K,  K-l)  system  would  be  the  value  of  K  such  that  K  =  p.i/p2  where  m  >  \i2,  and 
one  would  select  the  integer  closest  to  M-i/M-2  if  the  quotient  were  not  an  integer.  Us¬ 
ing  the  approximate  model  for  Q2  and  the  exact  results  for  Qi,  we  investigated  three 
cases  to  try  to  determine  if  this  were  so.  The  results  are  given  in  Figures  2,  3,  and  4. 
In  each  figure  the  values  of  E{Qi},  E{Q2)  and  E{Q}  =  E{Qi}  +  E{Q2)  are  given  as  a 
function  of  K.  For  all  three  figures  one  sees  that  the  minimum  of  E{Q}  does  not  oc¬ 
cur  exactly  at  K  =  p.i/p.2  but  is  close.  Figures  3  and  4  demonstrate  that  the  optimal 
value  of  K  does  depend  on  X.  as  well  as  pi/p2-  In  the  case  of  X  =  4  (pi/ji2=10),  the 
expected  number  of  customers  in  the  system  was  constant  for  K>5.  Whereas  when 
X  =  8  the  minimum  occurred  in  the  neighborhood  of  K  =  9. 


1  5  10  15  20  25 


FIGURE  2.  EXPECTED  QUEUE  LENGTHS  AS  A  FUNCTION  OF  K  F(JR 
1  X  =  3.5,  iJ.  AND  a2  =  1 

£ 
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.  3.  EXPECTED  QUEUE  LENGTHS  AS  A  FUNCTION  OF  K  FOR 
X  =  4,  AND 
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At  first  glance  one  would  expect  that  the  (K,  K-l)  scheme  would  result  in  a 
smaller  expected  number  of  customers  in  the  system  than  if  one  just  formed  one 
queue  and  placed  customers  in  service  as  a  server  became  free.  The  second  system 
was  studied  by  Singh  in  [2],  and  his  results  were  used  to  compare  with  the  multi¬ 
trunking  scheme  (K,  K-l).  This  analysis  is  given  in  Tables  2,  3,  and  4.  In  each  of  the 
tables,  the  total  service  rate  41  +  42  was  fixed  at  20  and  41  was  varied.  For  each  41, 
we  found  the  K*  which  minimized  the  expected  number  of  customers  in  the  multi¬ 
trunking  system,  and  we  used  the  EPP  approximation  at  system  2.  For  each  table,  one 
sees  that  the  expected  number  of  customers  in  the  multi-trunking  system  is  less  than 
that  under  a  single  queue  with  first-come,  first-served  rule  (Singh)  only  when  *11/42 
is  large.  But,  in  general,  the  results  of  the  two  rules  are  relatively  close  for  the  cases 
we  considered. 

As  we  pointed  out  in  Section  I,  the  multi-trunking  scheme  tends  to  keep  packets 
from  the  same  message  in  the  appropriate  order  of  arrival  at  the  receiving  node.  The 
main  reason  why  a  single  queue  with  first-come,  first-served  rule  (Singh)  sometimes 
results  in  a  shorter  total  queue  is  that  there  are  times  under  the  multi-trunking  scheme 
where  server  1  is  idle  and  server  2  has  customers  requiring  service. 

Several  other  interesting  facts  can  be  seen  from  these  tables.  First,  as  41  de¬ 
creases,  the  expected  number  of  customers  increases  for  the  multi-trunking  scheme. 
Thus,  if  41  ai  d  42  are  relatively  close  in  value,  then  one  would  not  want  to  consider  a 
multi-trunking  type  of  scheme.  For  Singh’s  system,  the  behavior  of  the  expected  num¬ 
ber  of  customers  in  the  system  is  convex  with  a  minimum  occurring  in  the  neighbor¬ 
hood  of  41  =  12. 
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M-l 

K 

Multi- 

Trunking 

19 

12 

f - " 

2.732 

3.129 

18 

8 

2.967 

2.988 

17 

6 

3.126 

2.888 

16 

5 

3.236 

2.818 

15 

5 

3.338 

2.771 

14 

4 

3.37  0 

13 

3 

3.481 

12 

3 

3.459 

11 

3 

3.507 

K 

Multi- 

Trunking 

Singh 

19 

19 

9.721 

9.847 

18 

14 

9.737 

9.678 

17 

11 

9.881 

9.636 

16 

9 

10.051 

9.567 

15 

8 

10.188 

9.517 

34 

7 

10.356 

9.484 

13 

6 

10.557 

9.464 

12 

6 

10.676 

9.457 

11 

5 

10.793 

9.460 

10 

5 

10.953 

9.474 

TABLE  4.  Expected  Numbers  of  Customers  -  Multi-Trunking  vs.  Singh  - 

A.  =  18  and  pi  ■*  42  =  20. 
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V.  Summary 

We  have  presented  exact  and  approximate  performance  models  for  a  multi-trunk¬ 
ing  scheme  for  packet  switched  networks.  The  exact  solution  was  numerically  compli¬ 
cated  by  the  fact  that  one  had  to  solve  for  roots  of  an  equation  and  then  solve  a  sys¬ 
tem  of  linear  equations.  The  approximation  required  none  of  these  complications  and 
gave  results  that  were  in  good  agreement  with  the  exact  approach. 

Using  the  approximate  model,  we  made  some  comparisons  with  a  single  queue 
with  first-come,  first-served  rule  and  found  that  the  multi-trunking  scheme  yields  a 
smaller  expected  number  of  customers  only  when  the  ratio  of  service  rates  is  large. 

Finally,  we  investigated  the  optimal  value  of  K  that  minimized  the  expected  num¬ 
ber  of  customers  in  the  (K,  K-l)  system.  It  was  found  to  occur  around  the  value  of 
K  =  ixi/m-2  (with  no  costs  considered).  A  very  interesting  problem  to  consider  would  be 
to  determine  whether  the  optimal  value  of  K  is  changed  under  an  imposed  cost  struc¬ 
ture.  For  instance,  if  there  is  a  significant  cost  to  switch  customer  arrivals  from  system 
1  to  system  2,  then  one  would  expect  the  optimal  value  of  K  to  vary  significantly. 


References 


[1]  M.  Steenstrup,  L.  Denenberg  and  H.  Mills,  “Store  and  Forward  Enhancement  Al¬ 
gorithms  Design,”  BBN  Report  No.  6360  (1987). 

[2]  V.  P.  Singh,  “Two-Server  Markovian  Queues  with  Balking:  Heterogeneous  vs  Ho¬ 
mogeneous  Server,”  Oper.  Res.  18,  145-159  (1970). 

[3]  S.  Halfin,  “The  Shortest  Queue  Problem,”  J.  Appl.  Prob.,  22,  865-878  (1985). 

[4]  C.  Knessl,  B.  J.  Mathowsky,  Z.  Schuss  and  C.  Tier,  “Two  Parrallel  Queues  with 
Dynamic  Routing,”  IEEE  Transactions  on  Comm.  VoI-Com-34,  No.  12  (1986). 

[5]  D.  Gross  and  C.  M.  Harris,  Fundamentals  of  Queueing  Theory ,  2nd  Edition,  John 
Wiley  &  Sons,  New  York  (1985). 

[6]  I.  L.  Mitrany  and  B.  Avi-Itzhak,  “A  Many-Server  Queue  with  Service  Interrup¬ 
tions,”  Oper.  Res.  16,  628-638  (1968). 

[7]  U.  N.  Bhat  and  M.  J.  Fischer,  “Multichannel  Queueing  Systems  with  Heterogenous 
Classes  of  Arrivals,”  Naval  Research  Logistics  Quarterly,  23,  271-282  (1976). 

[8]  B.  Kraimeche  and  M.  Schwartz,  “Analysis  of  Traffic  Access  Control  Strategies  in 
Integrated  Service  Networks,”  IEEE  Transactions  on  Comm.,  Vol-Com-30,  No.  10 
(1985). 

[9]  J.  W.  Cohen,  The  Single  Server  Queue.  Rev.  Edition,  North-Holland  Publishing  Co., 
Amsterdam,  (1982). 

[10]  A.  Kuczura,  “The  Interrupted  Poisson  Process  As  An  Overflow  Process,”  B.S.T.J., 
Vol.  52,  No.  3  (1973). 

[11]  H.  Heffes,  “Analysis  of  First-Come  First-Served  Queueing  System  with  Peaked 
Inputs,”  B.S.T.J.,  Vol.  52,  No.  7  (1973). 


-  24  - 


DISTRIBUTION  LIST 


